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Abstract 



We consider the quantum dynamics of a single particle in the plane under 
the influence of a constant perpendicular magnetic and a crossed electric poten- 
tial field. For a class of smooth and small potentials we construct a non-trivial 
invariant of motion. Do to so we proof that the Hamiltonian is unitarily equiva- 
lent to an effective Hamiltonian which commutes with the observable of kinetic 
energy. 



Consider a particle of mass m and charge e in the plane under the influence of 
a constant magnetic field of strength B and an electric potential. We choose 



the units of magnetic length y j^jj. the gyration time and the energy gap 

The dynamics are generated by 
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1 Introduction 




H = H La + V in L 2 (R 2 ) 



with 




with operator core the Schwartz space 5(M 2 ). V is the multiplication operator 
by a function V(q) and (qi,q2) ± '■= (— <Z2> (Zi)- For the gaussian 

g(q) := (q G M 2 ) 

we consider the class of functions defined by convolution with a real valued 
finite measure \i, g * p,(q) := f R2 g(q — q')dp,(q') : 

G := {V : R 2 -¥ R;V = g * ^, / d\fx\ < oo}. 

Jr 2 

Our main result is that there exists a non-trivial integral of motion; thus, in 
this weak sense, the two-dimensional system is integrable: 

Theorem 1.1 For V G Q small enough there exists a unitary operator U such 
that 

[U- l {H La + V)U,H La ] =0. 
In particular: UHloU^ 1 is an invariant of the flow e~ %m for all t G M. 

The meaning of "small enough" will be made precise in the sequel. A 
potential in Q is depicted in figure 1. 

It is folklore in plasma physics that in slowly varying fields the classical 
particle gyrates on a cycloid whose center drifts along the contour lines of the 
averaged potential and whose kinetic energy is an approximate invariant up to a 
certain time [CB, Nei]. On the other hand classically chaotic motion may occur 
in relevant regimes of parameters [GWNO]. 

In the realm of quantum physics the corresponding invariance is, in the 
large magnetic field limit, an essential ingredient for the current understanding 
of the integer quantum Hall effect [L, BESB, ASS, HS, GKS, CC, ABJ] while 
the interesting physics happen in the lowest Landau level. Coupling between 
Landau bands may, however, lead to non negligible effects [PG]. 

Several methods to construct an approximatively invariant subspace, making 
precise the notion of lowest Landau level, and an effective Hamiltonian are 
known, see [N, T, BDP]. They may lead to estimates valid to exponential order 
in a small parameter and valid for exponentially long times. See [BG] for an 
application of these ideas to a propagation problem. Our aim here is to study 
all Landau levels simultaneously as well as the possibility to go to the limit of 
all times. 

Our strategy is to employ a superconvergent iterative partial diagonaliza- 
tion procedure which was originally introduced in quantum problems to discuss 
stability of non-resonant time periodic problem [B, C, D$]; see [ADE] for an 
application to a condensed matter problem. It coincides in first order with 
the above mentioned "space adiabatic" algorithms. In particular the effective 
Hamiltonian restricted to the lowest Landau level is in first order (V)(x,D), 
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Figure 1: V(q) = J2ie[-io io] 2 u-i^diQ ~ *)> u (0 '-'-^ random variables 

a pseudodifferential operator whose symbol is the potential averaged over the 
Landau orbits. At higher orders the algorithm differs; on one hand it exhibits 
quadratic convergence, on the other hand it is an unsolved, and to our opinion 
important problem, whether our higher order effective operators are pseudodif- 
ferential. 

The partial diagonalization procedure is roughly described as follows: the 
diagonal part of an operator H with respect to a fixed, orthogonal, mutually 
disjoint family of projections {P n } is defined by DH := ^ n PnHP n . In fact, 
as the dimension of P n is not supposed to be one, or even to be finite, DH is 
block-diagonal. Suppose that the off-diagonal part OH := H — DH is small 
with respect to DH. The equation 

e w He~ w -DH = 0(\\OH\\ 2 ) 

is satisfied by an antiselfadjoint operator W of order O (||OiT||) solving 

OH + [W, DH] = 0. 

If one takes P n the projections on the Landau levels then, because of the gap, a 
solution of this equation can be found if the coupling between the bands decays 
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sufficiently fast. This is the case for the potentials of our class Q. The procedure 
can then be iterated by replacing H by e w He~ w . The convergence of the 
transformed H to a block diagonal operator which, because of the degeneracy, 
commutes with PiLa is quadratic. 

We remark that in order to be really applicable to the quantum Hall effect 
our result should be extended to a more general class of potentials than stated 
above. This should in principle be possible. However, a delicate control of 
[HL a ,V] is needed. The method does not work for the extensively studied 
purely periodic problem V(q\,q2) = cosgi + cos (72- 

The plan of the paper is to set up the iterative algorithm in section 2. The 
class of potentials Q and control the necessary norms will be discussed in section 
3. Section 4 contains the proof of theorem 1.1 and a discussion related special 
cases. 



2 An iterative partial diagonalization algo- 



As discussed in the introduction, the task is to partially diagonalize the operator 



Recall that Hi a = ^neN ( n + V^) w ' tn infinite dimensional projec- 
tions P^ a ■ We consider H which is of the same type as H^a, in order to cover 
situations where H is already an effective Hamiltonian at finite order we assume 
that it has either a finite number of bands a n or an infinite number such that 



Definition 2.1 A selfadjoint operator H is of class C s for a g > if for a 
complete family of orthogonal, mutually disjoint projections {P n }neIcN which 
commute with H, it holds for 



rithm 



H La + V. 



dist (a n , a m ) > g|n — m 




n 



RanP, 



n 




min<7 n +i — max a n > 0- 



... [. 



E — ] 



1- 



Figure 2: Typical spectrum of H 
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For the same family of projections P n and a bounded operator V we use the 
notation 

DV := PnVP n , OV := V - DV. 

To organize the estimates we shall make frequently use of the notations : 

(a) := max(l, |a|), || || : the operator norm, and 

II^Hj := sup {n - m) l \\P n AP m \\. 
n,mEl 

We prove (extending [DLSV]) 
Theorem 2.2 Let H G C s and V be a bounded selfadjoint operator such that 

mu < |. 

Then there exists a unitary U such that U^ 1 D(H) C D(H) with the property 
that for 

H^:=U{H + V)U-\ D(H 00 ) = D(H) 

it holds 

[Hoo,P n ] = 0. 

Proof. Define H := H + V = UH + OH . Assume ||V||i < |. Then 
DHo G C1/4, [|Oflo[|i = ||OV[|i < oo. By lemma 2.3 there exists a bounded 
solution Wo of 

[T>H ,W ] = OH , dw = o. 

Define Wo := e Wo . 

Remark : U is unitary, D(H ) C D(H), W D(H ) C D(H ), U D(H ) C 
D(H ). 

Suppose that for s G N diagonalization has been done up to H s , U s -\ such 
that ||Di? s — H\\ < 0/4, ||Oi? s ||i < 00. To go to step s + 1, use lemma 2.3 
to solve 

[Dtf s ,W s ] = Otf s , W s = 
for a bounded W s and define U s := e^Z/g-i, 

tf s+1 := e w 'H.e- w ' = e L ^(H s ) = ]T 

fc=0 

with the notation L W (A) := [W,>1]. Now 

%(Dff s ) = -OH s 

thus for A; > 1 

L k Ws (H s ) = -L k wHOH s ) + L k Ws {OH s ) 
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so 

H s+1 = -DH s + 4>(L Ws )(OH s ) (1) 

with 

1 °° k 

4>(x) := - -{? - l) = E ^7^* (* e R). 
It holds by corollary 2.5 

11x^(0^)111 < ||w.ofl a ||i + ||Oiw s ||i < 2^11^112110^111. 

Thus by induction and lemma 2.3 

k 

k+1 



( tt2K\ 

\LUOH s)h <^—) \\OH s \\\ 



which implies 



\\<KLw.)(OH B )\\i < E(fcTl)! (-g-H^lli) l' 0// - 
- ^||Ofl s ||i) 



2vrA" r v g 

with tp(x) := x<ft(x) = (x— l)e x + l (x > 0). Remark that -0 is nonnegative, 
< x (x < 1), and that is a superattractive fixed point. By (1) 

Off^ i < ^ I Ofl s | h 



thus 

l|Ofl s ||i < A 2 * (2) 
with A proportional to ||OV[|i small enough. For the diagonal part it holds 

a 2-kK 

HDflH.i-DH.il = ||Dfl s+ i-Dfl s ||i < ||(^(X w J(Ofl s )||i < l|Ofl s ||i) 



and thus for an x proportional to ||OV|| 



||Dfl s+ i - fl|| < * + ^< Q - 
11 + 11 - 8 2-kK J y ~ 4 

so the iteration is well defined. 

In particular s h-> Dfl s — H is a norm convergent sequence. (2) and Lemma 
2.3 imply that ||W S || converges superexponentially to 0. By (1) this implies in 
turn as OV is bounded and ||OV||i is small enough: 

||Ofl s+ i|| < H(L Ws )(OH s )\\ < E J]^r[y2 k \\Ws\\ k \\OH.\\ 

V v ' 

<1 
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thus OH s converges in operator norm and furthermore as 
||W s +i — W a || = ||e M/s+1 — I|| one concludes that U s — U. By construction 
= D^/qo which commutes with the projections. H 
We now prove some claims which where used in the preceding proof: the 
existence of a bounded solution of the commutator equation is assured by 

Lemma 2.3 Let H £ C g and V a bounded selfadjoint operator such that 
PnVPn = Vn and such that \\V\\\ < oo. Then there exists a bounded 
antiselfadjoint W such that 

[H, W] = V, BW = 

such that 



\\Wh < £-\\V\\i. 

Proof. By [BR] for bounded operators A,B,C there exists a solution X of 
AX - XB = C such that 

Ul < * 11011 



2 dist (spect(A), spect(B)) 
It follows that there exists W nm such that 

PnHPjivV nrn W nm P m H P m — P n V P m 

with W nm = P n W nm P m and 

7T ,. 1 7T ||V||l 

\\W nm \\<-\\P n VP m \\-, r< " " 



2 n m g(n — m) 2g (n — m) 2 ' 
Now define W := J2 n ^m,Wnm in norm convergence, then 

m s eupEiiw».ii - ^E („_' m)1 "'"' 

from which the claim follows. H 

Lemma 2.4 Z.et n,m G N anc/ K := 3 + 2£(2). /t /jo/cfs: 

1 < ^ 

■^n^-^ {j ~ n) 2 (j - m) ~ {m- n)' 
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Proof : The case where n = m is evident. In what follows, we will simply 
write j as the index for the sum except of j > 0, j ^ n,j ^ m. 

1 _ 1 (\r^ {j-n + m-j) 

{j ~ n ) 2 {j — m ) ( m — n ) \~ J — n ) 2 {j ~ m ) 

1 1 v 1 

~~ (m — n) \ (j — n) 2 ^ {j — n) {j — m) 

The left term of the r.h.s. is bounded by 2£(2). It remains to prove that 
the second term of the r.h.s. is bounded. Recall the series expansion of the 
digamma function Vo : 

oo 

V, (x + l) + 7 = ^— ^— ( xgN )- ( 4 ) 
~{JU + x ) 

Now, since it is symmetric in m and re, we can assume that m > n and define 
a := m — n > 1. Then the second term of the r.h.s. becomes : 

y- L = y 1 



, (j - n) (j - m) .j^ n (j) (j - a) 



a-l 



V 1 v 1 

(j')O'-a) + (j)(j~a) 

j=-n w/w ' 3>a+l w/w ' 

1 7 + ^o(l + a) 

ji^tiW-*) a 

where we used (4). But 

a— 1 ^ n ^ a— 1 ^ 

^Jj)(j-a) = ]^(j)(j + a) + ]^(j)(j-a) 

7 + ^o(l + a) 1 

< Hfl-1 SUp —r- r. 

a je[i,a-i]nN - a| 

Since a > 1, both of these terms are bounded by 1 so the claim of the lemma 
follows. I 

We have the following corollary : 

Corollary 2.5 For operators A and B such that \\A\\2 < oo and \\B\\i < oo it 
holds 

\\AB\U < KWAMBh 
where K was defined in 2.4. 
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Proof : Since {P n } n i=i is a complete family of orthogonal and mutually 
disjoint projectors, we can write 

P n ABP m = Y J (PnAP l )(P l B 

Pm)- 

lei 

Thus 

\\AB\\ 1 < S u l (rn-n)^^ 2 A 
The result follows from lemma 2.4. I 



3 The class Q 

We show that the basic decay estimate is satisfied for potentials in the class Q 
defined in the introduction and give some examples. 

Proposition 3.1 For V G Q and P n the eigen projections of Pha on its n-th 
level it holds in operator norm on L 2 (M 2 ) for a d > and all n,m 6 No-' 

II p VP II < 



(n — m) 



Proof. 

\\P n vp m \\ = \\PnJ 9 (.-y)My)P m \\ 

< fw(y)\ \\P n g(. - y)P m \\. 

Consider the unitary magnetic translations on L 2 (M 2 ) defined for a a G M? by 

T{a)t/j(q) = e^ gAa ^(q-a). 

It holds [T(a),H La ] = and T(a)gT*(a) = g(. - a). Thus 

\\P n g(.-y)P m \\ = \\P n gP m \\ Vy € M 2 

and the result follows from Proposition 3.5 to be proven below. I 

Remark that V = g*fi extends necessarily to an entire analytic function; its 
Fourier Transform V has the property that V{p) exp \ is the Fourier Transform 
of a finite measure. We elaborate on this in order to point out that Q contains 
sufficiently many potentials to be of interest for applications to the Quantum 
Hall effect. 

Definition 3.2 Denote by A the class of functions V : M. 2 — > R such that 
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1. V has an extension to an entire function on C 2 , 

2. R 2 3 y h-> e-Vv(iy) G L 1 ^ 2 ), 

3. for K(p) := e 2 ? / R2 e'^e'^ V{iy)j0p it holds: V G L^R 2 ). 
Proposition 3.3 For V £ A it holds 

V = g*(Vdq). 
Proof. By Fourier's theorem, it holds for gel 2 



,2 



dy 



e -^V(iq) = I e v(5-v) e -Vy(iy)-^-dp. 

ixR2 I 2 ?!") 



Thus 



V(it/) = / e 2 V(p)dp 
Jm. 2 

and the claims follows as both sides are analytic in q. H 

We list some examples of potentials in Q which contain in particular Ander- 
son type models on a finite portion of the probe. 

Corollary 3.4 Q contains 

1. for p a real polynomial, a G (0, 1), k\, ki G R: 

p{qi,q2,e ikiqi ,e ik2q2 )e- aa 2 2 -, 

2- E ie z 2 Vi9{q - i) with ii G l\l?,R). 
The function q H> e lkq for k G R 2 does not belong to Q. 

Proof. For 1. it is sufficient to verify that 

y 2 i 
R 3 y i-)- e 2 f(i y ) e L 

and 



2 r , 2 
Bn-Je^ / e~ lxy f (iy)dy G L x ( 



for /(y) = y n e~ a V and f(y) = e iky e~ aM 2 , „ G N and fc G R. 

This is by standard properties of the Fourier transform. In the first case, 



\f(iy)e- ay -T\ = \y\ n e-^- a ^ eLH 
and 

f e- ixy y n e- (1 - a) £dy\ = . 1 \poly{x)e^ e~^^ | G L 1 (R). 
J Vl - a 
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si 
e 2 



In the second case 

\f{iy)e~ a £ \ < e \ kM e- {1 - a)V ^ € L 1 (M) 

J- f i n w 2 1 fc 2 x 2 1 x 2 . 

\e* / e- lxy e- ky e- {l - a) - 2 dy\ = p W^) p ^ p ~t^^ G L 1 ^) 

y Vl — ol 

In 2., one deals with the pure point measure 

Vi5{x - i). 

Finally, one has 

e ikx = g*v 

with /i = e e 2 <ix which is not a finite measure. ■ 

It remains to treat the case of the gaussian potential g which turns out to 
be non trivial. We follow the strategy designed in [W]. 



Proposition 3.5 For g(q) = exp ( — ^-J, (q G R ) and P n t/?e eigenprojec- 

tions of H^a on its n-th level it holds in operator norm on L 2 (M 2 ) for a c > 
and all n,m G No-' 

-fn^-^m — "7 r- 
(n — m) 

Proof. We use the representation of P n by angular momentum eigenfunctions 

P n= ^2 \^ri,l){^n,l\ 
l>-n 

1 1 / 2X - r - 



Mr,S);=(-irJ ¥ ^r' e «L^-)- s (5) 



where the Laguerre polynomials are defined by 

:= (^(-x)WlJ^(x) (0>I>-n). 



Then 



thus 



P n gPm= ^ I^W)(VV'0Vw)(VV, 



l>— nAm 



|(t/',Pnfl r -Fm¥')| < sup | (V> n ,z , WW 

i>-nAm 

and the claim follows from equation (6) and estimate (9) to be proven in the 
following two lemmas. H 
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Lemma 3.6 For g(q) = exp y— \ J , (q G R 2 ) and ip n j defined in (5) it holds 

for n,m G No, I > — (n A to) 



i// / \i 1 (Z + m + n)! 

ft,i.#m,l = „ f+m+n+1 -7 ., ■ .... ■ Ff=j , - (6) 



Proof. By definition 



Consider first / > 0. To study the dependence of the integral on l,m,n we use 
that the family n \— > L n (x) is orthogonal in L 2 (S. + ,di , i), dv\ := x l e~ x dx and 
the identity: 

^n^m (2) = ^ ' m ' 1 L\ (x) . 
s>0 

As L l = 1 VZ, one has 



poo r poo 

/ L l n l) m {-)du l {x) = Y,B^ 1 L\l^{x)dm{? 
Jo * s > Jo 

poo 

= B% m ' 1 dv l {x) = B^ m ' l T(l + l) 
Jo 

It was proven in [Ca] that 



9s(x,y,l) 



1 v_ 



2(l-x) ^ 2(l-j/) 



(l-x)'+l(l-y)I+l (l + ^ + ^) 



/+a+l 



Thus 



Bl ' m ' 1 = J^a?qgo(x, y, I) 
rurrv. y 



1 (l + m + n)\ 
=0 ~ 2 m+n l\m\n\ 



from which the claim follows for I > 0. 
Now for I < 0, one has 

VV = ^n+l-l 

thus 

\(lpn,l,Vlpm,l) \ = \{lpn+l,-l,Vlp m+ l-i)\ 

and the result follows for I < 0. I 

The following lemma might be known to probabilists, we know of no refer- 
ence though. 
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Lemma 3.7 For a c > 0, it holds 

(m + n)\ m + n ™- n , , . 

< c =e 2 < m +") (N9m,n>l), (7) 



2 m+n n!m! - + n )2 _ ( m _ „)2 



(m + n)\ c 

<-, r (m,ne N ), (8) 



2 m+n n\m\ (m — n) 

(1 + m + n)\ c . , 

< 7 r (m,neNo,l>-mAn). 



gz+m+n+i v / (/ + m)!(/ + n)!n!m! (m - n) 

(9) 

Proof. Denote by C a positive constant whose value may change from line to 
line. 

We use Stirling's and a concavity inequality : 



1 < j^t^ <C (n > 1) (10) 



(i 

(1 + x) log(l + x) + (1 - x) log(l - x) > x 2 (arG [0,1]) (11) 

to estimate 

(m + n)! 
flm ' n := 2™+™n!mf 
For m,n > 1 it holds by (10) : 



( m + n )(m+n) /J I 



<hn.n < C — - — zn^—rz. \ — I < C 



m ' n - m m n n 2 m +«Vn m~ ( a + &\*±* ( a _ V a 2 - 6 2 

with a := m + n, b := m — n in 

G := {(a, 6) G Z 2 ;a > 2, |6| < |o-2|} . 

Note that the case n = m follows from the first inequality of the previous line. 
Using (11) with x := ^ it follows 



<l,n,n < \ I q2 _ ^ 2 e -■' 

which implies (7). Consider now 

[ a z 2 I 

then 2 

(m-n) 2 a 2 rti „ < C-^~e'£ < 2 ((a, 6) € G<). 
1 — -75- a 
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In G\G<, it holds 

a b 2 n b 2 2 

< — < a as well as 1 > - 

2 a a z a 

thus 

(m — n) 2 a 2 mn < C a 2 e~ 2 < C 

which proves (8) as it is evident for n = 0. 
Denote now 



bl m,n •" 



(l + m + n)l\ 2 1 



2 l+m+n J (J + m )\(l + n )\ m \nV 

From the identity b^ m ^ n = ai +rriyn ai +n ^ m it follows for I + m, I + n, m, n > 1 



with a := I + m + n,b := I + m — n, c := I + n — m, (a,b) G G and (a, c) G G. 
It follows 



(m - n) 2 bi tmtn < C 



b-c 



>a 



2 



1 



b 2 



-e 2a . 



In G, one has — ^ < | so if ^ > | then | < ^ < a and, denoting any 

1 2 

polynomial by "po/y" : 

C — a c 2 

(m - n) 2 b hm ^ n < Cpoly{—=, y/a)e~%e~^ < C. 
Now if < i then either 4" > ij and 

or 2 ar — 2 

_ b 2 +c 2 

(m - n) 2 bi^ n < Cpoly(^, -p)e ^~ < C 

V a 

or ^ < 1 and 

& C b 2 | 

(m - n)\ mn < Cpoly{—=, —=)e ^ < C. 
For the cases where one of I + m, I + n, m, n is zero, remark first that 

2 c 

bl,n,n = | a Z+n,n| < "^p" < C 

by (8) and secondly that for I + m = 0, n > m 

2, 1 ,9 (m — n) 2 C 
(m — n) bi m „ = — (m — n) a/+„ m < ; r < C 

which covers all cases. H 
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4 Application of the algorithm 



We proof Theorem 1.1 then give an illustration. 
Proof, (of Theorem 1.1) 

Choose H = Hi a and {P n } nl =j$ its eigenprojections. Then Hi, a G C\. By 
proposition 3.1 ||V[|i is finite for V G Q. So for a V G Q with < \ by 

theorem 2.2 there exists £Y unitary such that [U (H^a + V) U~ l , P n ) = thus 
[U(H La + V)U-\H La }=0. I 

4.1 Quadratic Hamiltonians 

We discuss the case where V is a polynomial of degree at most 2 for sufficiently 
high magnetic field. Though this case is not covered directly by Theorem 1.1, 
the iterative algorithm can be applied to the hamiltonian matrix which defines 
the operator. This results in the construction of an integral of motion which is 
the quantization of a classical integral, independent of the hamiltonian function. 
The following operations are to be understood first on vectors in 5(M 2 ) then 
on the appropriate extensions. Denote D := — iV, the velocity and center 
operators 



and recall the commutation relations 

[v!,v 2 ]=i, [c 2 ,ci]=i, [ci,Vj}=0. 
The linear case is trivial, nevertheless it is instructive : 

V(q) = -(E,q) = -(E iqi + E 2 q 2 ), 
define W := i(E,v), Uq = e w °. From the Weyl relations 

e i(E,v) ve -i(E,v) =V _ E ± 

it follows 

Proposition 4.1 For E G R 2 , V{q) = ~{E,q), U = e*< E ^ , it holds 

1. U (H La + V)Uq 1 = H La - (E, c) - \E 2 = BH- \E 2 

2. U Q 1 H La Uo = U v + E± ) 2 - 

Now consider the quadratic case, V(q) = \{q,V"q) for a real symmetric 2x2 
matrix V" . 

The Hamiltonian is 
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with 



a t /2 -I 
1/4 + V" a*/2 



1 N 

where we denote a = ( ] and I the 2x2 identity matrix. 



is a hamiltonian matrix with respect to the symplectic structure defined 

I 0, 



, i I i i 

by . Its eigenvalues are 



{±iy/P + Q,±y/Q-P} 



with P := 1 + trV", Q = P 2 -4det V". For V" small enough, P > 0, Q > 
thus A G zM. In the case of an hyperbolic fixed point, det V" < so /x € M. 
In the elliptic case, /j G E if det V" is small enough and in the parabolic case, 
/i = 0. In all three cases one knows from normal form theory (see [Wi]) that 
there exists a symplectic transformation decoupling the degrees of freedom. We 
state the explicit result for the cases of the quantum dot and antidot, i.e.: 
V" = ±e 2 I which one verifies by direct calculation: 

Proposition 4.2 Let 1/2 > e > 0, V{q) = ±^(q 2 + q%). Define ft := 
Vl ± 4e 2 and the unitary Uip(q) = -^ij)(q/y/Q), then it holds 

1. U(H La + V)U~ X = ^H La + «=lf = Dff + ^(H La + f) T 

2. anc/ /"or t/?e constant of motion: 

i 
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